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We describe a simple model for the melting and optical properties of a DN A/gold nanoparticle aggregate. 
The aggregate is modeled as a cluster of gold nanoparticles on a periodic lattice connected by DNA bonds, and 
the extinction coefficient is computed using the discrete dipole approximation. The optical properties at fixed 
wavelength change dramatically at the melting transition, which is found to be higher and narrower in temperature 
for larger particles, and much sharper than that of an isolated DNA link. All these features are in agreement with 
available experiments. 



1. Introduction 

Recently, so-called functional metallic nanopar- 
ticles have started to be developed, which may 
lead to new materials with improved optical and 
mechanical properties [1] . Among these, there is a 
particular interest in gold nanoparticles to which 
noncomplementary oligonucleotides capped with 
thiol groups are attached (DNA modified gold 
nanoparticle system), because in addition to a 
strategy using self-assembly of nanoparticles [2] , 
biological detection is possible using the optical 
and electrical sensitivity of their aggregates [3,4]. 

In this paper, we present a model for this "melt- 
ing" transition, which accounts for most experi- 
mentally observed features. 

2. Structural Modeling 

To calculate the T-dependent optical proper- 
ties, we have considered two slightly different 
structural models. In the simplest version of our 
model, the low-T aggregate is taken simply as 
a collection of identical gold nanospheres (each 
of radius a) which form a cube of edge L and 



fill the sites of a simple cubic lattice of lattice 
constant d (d > 2a), containing iVpar = {L/d)^ 
gold nanoparticles. We have also investigated the 
melting and the optical properties assuming that 
the low-T cluster is a fractal aggregate. In both 
cases, we assume that all bonds are occupied by 
the same number [Ng/ z) of DNA links at low 
temperature, where Ng is the number of single 
DNA strands attached to one particle and z is 
the number of nearest neighbors (z = 6 for a sim- 
ple cubic lattice). 

As temperature is increased, the bonds start to 
break. We define the probability that a link forms 
between a given gold particle pair at temperature 
T as Pcr{T). To generate a specific sample at 
temperature T, we randomly remove links with 
probability 1 — pcff(r), then identify the separate 
clusters, using a simple computer algorithm [5]. If 
the aggregate is separated into two or more clus- 
ters, we simply place these clusters in random po- 
sitions and orientations within a larger bounding 
box (usually of edge lOOd), taking care that the 
individual clusters do not overlap. The resulting 
geometry is shown schematically in Fig. 1. As 
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a second structural model, we have carried out 
the same procedure to simulate the melting of a 
sample formed by reaction-limited cluster-cluster 
aggregation (RLCA) [6]. A typical RLCA clus- 
ter is shown in Fig. 1 (d), and represents a pos- 
sible fractal aggregate which might be produced 
by certain random growth processes at low tem- 
perature [7]. Here we implicitly assume that the 
specific bonds which are occupied at temperature 
T are time-independent. 

3. Melting Theory for DNA Links 

The two-state model can describe melting well 
for short DNA strands (12-14 base pairs) [8,9], 
and is given by the relation 

S + S^D. (1) 

The chemical equilibrium condition correspond- 
ing to (1) is 

[1-P{T)? _ K(T) 
P{T) - Ct ' 

where p{T) is a (tcmpcraturc-dcpcndcnt) frac- 
tion of the double DNA strands, K(T) is a suit- 
able chemical equilibrium constant, and Ct is 
the molar concentration of single DNA strands 
in the sample. We have also assumed the sim- 
ple van't Hoff behavior K{T) = exp[-AG/kBT], 
with a Gibbs free energy of formation AG(T) = 
Ci{T - Tm) + C2{T - Tm)^ , choosing the values of 
ci, C2, and Tm to be consistent with experiments 
on these DNA molecules. 

In order to calculate the fraction Pes{T) of 
bonds which contain at least one double strand, 
we adopt the following very simplified model. 
First, we assume that exactly Ng/z of the sin- 
gle strands on a given nanoparticle are available 
to bond with any one of its z nearest neighbors. 
Thus, we assume that the maximum number of 
links that could be formed between any two par- 
ticles is Ns/z. The probability that no link is 
formed is then taken to be 

l-Peff(T) = [l-p(^r•'/^ (3) 

In Fig. 2, we plot Peff(T) for several radius a. 
We assume Ng oc a^, set z = 6 and use the exper- 
imental result that Ng = 160 when a = 8 nm [10]. 



In actuality, there is a linker molecule which 
emerges from solution to connect two DNA sin- 
gle strands on different nanoparticles. But in fact, 
we can avoid considering the linker molecule ex- 
plicitly, as shown in the inset of Fig. 2. In this 
inset, the dotted and dashed lines represent p{T) 
for two melting curves, one for each single strand 
connected to a linker molecule with slightly dif- 
ferent constants Ci and C2 but the same melting 
temperature. (This is the case where the linker 
molecule is expected to have the most effect.) The 
solid line is the rcsidting melting curve for the 
two single strands plus linker. The squares are 
the melting probability p'{T) calculated assum- 
ing two effective strands and no linker strand with 
slightly different ci, C2, and Tm- 

4. Optical Properties 

We calculate the optical properties of this 
sample using the Discrete Dipole Approximation 
(DDA) [11,12]. The sample is modeled as a 
collection of separate aggregates, whose extinc- 
tion coefficients are computed individually, then 
added. Each aggregate consists of many identi- 
cal nanoparticles, which have complex frequency- 
dependent dielectric constant e{Lo), and polar- 
izability 0(0;) related to e{u)) by the Clausius- 
Mossotti equation with radiative reaction correc- 
tion term [13]. The resulting expressions for the 
induced dipole moment Pi of the i*^ sphere, and 
the corresponding expression for the extinction 
coefficient Cext{k) at wave number k, are given in 
Ref. [14]. 

In our case, each cluster consists of a number of 
DNA-linked individual gold nanoparticles. In our 
calculations, we do not include the optical prop- 
erties of the DNA molecules, since these absorb 
primarily in the ultraviolet [15]. We use tabu- 
lated values of the gold complex index of refrac- 
tion [16,17], then calculate Cext for each cluster 
using the DDA. To improve the statistics, we av- 
erage Cext for each cluster over possible orienta- 
tions. We then sum the averaged extinction coef- 
ficients of all the individual clusters to get the to- 
tal extinction coefficient of the suspension. This 
method is justified when the suspension is dilute. 

In Fig. 3, we show Cext (A, T) at fixed A = 520 
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nm, close to the isolated-particle surface-plasmon 
resonance (SPR), versus T, for several particle 
sizes, assuming that the low-T aggregate is a sim- 
ple cubic cluster with TV = 1000 particles, which 
is illustrated in Fig. 1 (a). As can be seen in 
Fig. 3, for each radius a, the extinction increases 
sharply at a characteristic temperature, corre- 
sponding to the melting of the aggregate for that 
size; at this temperature, the absorption due to 
the SPR increases sharply. 

In the inset of Fig. 3, wc compare Ccxt (A, T) for 
a regular and an RLCA cluster of particles of 20 
nm radius at A = 520 nm, both for N = 1000 
particles. Although the RLCA cluster has a 
slightly broader melting transition, as manifested 
in Cext(A,T), than does the regular lattice, both 
sets of data show a much sharper melting tran- 
sition than that of a single DNA link. Also, al- 
though our normalized Cext(A, T) is calculated for 
the aggregates at 520 nm, we expect similar be- 
havior at 260 nm. (Wc have not carried out cal- 
culations at this A mainly because we have not 
included the DNA absorption properties.) In any 
case, the experimental melting curves at 260 and 
520 nm are very similar [18]. 

5. Discussion 

We have briefly described a model for struc- 
tural development of DN A/gold nanoparticle 
clusters, and have calculated the resulting clus- 
ter optical properties using the DDA. 

Our calculated extinction coefficients are good 
agreement with recent experimental results [3,15, 
18]. In particular, both experiment and calcula- 
tion give a sharp increase in Ccxi(A,T) at fixed 
A, as T increases past a critical temperature. 
We also find, in agreement with experiment, that 
melting occurs over a much narrower range of T 
in the aggregate than for a single bond, and that 
the melting occurs at higher T for larger particles. 
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Figure 1 . Schematic of the melting of a gold-DNA 
cluster, for two different models discussed in the 
present paper. In the first model, (a) shows the 
cluster at low T (pcs = 1); (b) peff = 0.5 > Pc{L), 
where Pc is the weakly L-dependent percolation 
threshold; (c) Pes{T) = 0.2 < Pc{L); (d) alternate 
model for low-T sample {pctf{T) = 1): fractal 
cluster formed by reaction-limited cluster-cluster 
aggregation [RLCA], with fractal dimension rf/ ~ 
2.1. 



carried out using the facilities of the Ohio Super- 
computer Center. 
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